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ABSTRACT 

We investigate effects of the magnetic field on the propagation of laminar flames of nuclear 
reactions taking place in white dwarfs (WDs) with the mass close to the Chandrasekhar limit. 
We calculate the velocities of laminar flames parallel and perpendicular to uniform magnetic 
fields as eigenvalues of steady solutions for magnetic hydrodynamical equations. As a result, we 
find that even when the magnetic pressure does not dominate the entire pressure it is possible 
for the magnetic field to suppress the flame propagation through the thermal conduction. Above 
the critical magnetic field, the flame velocity decreases with increasing magnetic field strength 
as v ~ B . In media with densities of 10 7 , 10 8 
orders of ~ 10 10 , 10 , andlO 12 G, respectively. 



and 10 9 gem , the critical magnetic fields are 
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supernovae: general — white dwarfs 



Type la supernovae (SNe la) are thought to be 
thermonuclear explosions of C +O white dwarfs 
(WD s ) in close binary sy stems (jlben fc Tutukov 



1984 iNomoto et al.|[l983) . If the mass accretion 
rate from the companion star is ~ 10 _7 iWQyr _1 , 
the WD can increase its mass to the Chan- 
drasekhar limit (1.38 M@) and the evolution ends 
up with the ignition of carbon in the innermost 
region. As a result, the explosion is expected to 
be regulated by this limit for all of this type of 
SNe to behave in a similar manner (ideally, in an 
identical manner). Thus they are used as standar d 
candles (|Riess et alJll998t|Perlmutter et al.lll999l) . 
In reality, observations of SNe la suggest that their 
peak luminosities h ave a significant diversity (e.g. 



Hamuv et al. 19961) . Despite of this diversity, we 



can estimate the absolute magnitude of maximum 
light from the decline rate of the luminosity after 
the peak using a tight correlation between them 
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(|Philhpsl 119931) . With this calibration method, 
SNe la are used as a cosmological distance indi- 
cator, which r eveals the accelerating expansion of 



the U niverse (jRiess et al 
1999h . 



1998; Perlmutter et al 



Yet, we do not agree on the intrinsic mech- 
anism for S Ne la to generate t h eir lu minosity 
dispersion. Hoeflich fc Khokhlov ( 19961 ) investi- 
gated different explosion scenarios including de- 
flagrations and detonation models and compared 
their light curves with those of observed SNe la. 
They found a correlation between the types of 
host galaxies and explosion mechanisms. The 
metallicity in a progenitor star is suggested to 
affect the composition of the WD and influence 



the b r ightness of the resul ti ng SN i Hocflich et al 



19981: lUmeda et all Il999t iTimmes et all 120031) 



The double degenerate scenario, in which a bi- 
nary system composed of double WD merges to a 
massive WD with the mass exceeding the Chan- 
drasekhar limit to result in an explosion, naturally 
yields a v ariety of explosions dep ending on the to- 
tal mass (jlben fc Tutukovlll984l ). Another possi- 
bility is the effect of mag netic field in a progenitor 
WD (jGhezzi et all 12004 ) ■ The propagation of ex- 
panding flow generated by deflagration may be 
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suppressed in the direction perpendicular to the 
magnetic field, which will cause an asymmetric 
explosion. Thus a diversity of SNe la caused by 
different strengths of magnetic fields depends not 
only on the amount of synthesized 56 Ni but also 
on the lines of sight of observers. 

Recent observations, especially Sloan Digital 
Sky Survey (SDSS) have discovered a number of 
WDs with strong surface magnetic fields greater 
thanlQ 4 G (ISchmidt et alJ2003 : IVanlandingham et al 
2005). A significant fraction of them have polar 
magnetic field strengths ranging from 10 6 G to 
10 9 G. There is a tendency that more massive 
WDs have stronger magnetic fields. Furthermore, 
WDs in close binary systems tend to have strong 
fields, which indicates that suc h magnetism is pro - 



duced in common envelopes (|Tout et al.l [2008) 



Note that these field strengths are the surface val- 
ues. At the center of a WD, it would be stronger. 
Since the mass of observed magnetic WDs are sig- 
nificantly less than the Chandrasekhar limit, (the 
average mass is ~ 0.9 M ), increasing their masses 
toward the Chandrasekhar limit, their fields are 
likely to become much stronger when these mag- 
netic WDs explode as SNe la. 



Ghezzi et al.l (|2004l ) studied effects of the mag- 



netic field on the turbulent combustion front and 
calculated resulting asymmetry of the explosion by 
evaluating the growth rate of the Rayleigh- Taylor 
instability of the turbulent front suppressed by the 
existence of magnetic fields. 

The effect of magnetic field is not only dy- 
namical. Electrons are trapped in the presence 
of the magnetic field, which suppresses the trans- 
verse thermal conduction. The velocity of a lam- 
inar defl agration depends on th e thermal con- 
duction (jLandau fc Lifshi"t3 19591 ). In the react- 
ing region of a deflagration front, the diffusion 
timescale Tdis = 1 2 /k, and the burning timcscalc 
Tbum = E/S arc comparable. Here, I is the mean 
free path of electrons, n the diffusion coefficient, 
E the energy density, and S is the energy gen- 
eration rate. Equating these timescales, the flame 
velocity relative to the unburned medium is ap- 
proximately estimated by the following formula, 



(1) 



This suggests that the transverse laminar flame 




slows down in the presence of the magnetic field, 
because the diffusion coefficient is proportional 
to the thermal conductivity. Thus, this effect is 
needed to be taken into account for discussing pos- 
sible mechanisms to give rise to an asymmetric ex- 
plosion and the diversity of SNe la due to magnetic 
fields. 

In this paper, we study the propagations of lam- 
inar flames in magnetic C+O WDs. Since the 
width of the wave front is much smaller than the 
scale lengths of a WD, we consider flames propa- 
gating into uniform media with constant densities 
and uniform magnetic fields. In this approxima- 
tion, we evaluate the critical magnetic field to start 
affecting the flame velocity for each density. 

The structure of the paper is as follows. In 
Section 2, we present formulation for the steady 
laminar nuclear flame. Section 3 explains the nu- 
merical procedure to solve the equations satisfying 
boundary conditions and to obtain the propaga- 
tion velocity as an eigenvalue problem. Section 
4 presents results of the calculations. Section 5 
concludes this paper. 

2. Formulation 



The laminar flame velocities in non-magnetic 
WDs are calculated bv lTimmes fc Wooslevl (|1992t ) 
(TW92). They integrate time dependent one- 
dimensional hydrodynamical equations for lami- 
nar flames propagating in uniform media till they 
obtain approximately steady states. Then the ve- 
locity of the front is measured. 

Instead, we integrate time independent one- 
dimensional steady state equations in the wave 
rest frame with a mass flux chosen so as to satisfy 
boundary conditions imposed at both ends of the 
flame front. This method allows us to determine 
the velocity of the wave front as an eigenvalue of 
the steady solution. 



2.1. MHD equations 

Since we consider one-dimensional steady flows, 
physical variables including magnetic fields de- 
pend only on the spatial coordinate x. MHD equa- 
tions governing the flow perpendicular to the mag- 
netic field can be written in the following manner. 
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where T is the temperature, p the density, v the 
velocity, P the pressure, h the specific enthalpy, 5" 
the net energy generation rate including the heat- 
ing process by nuclear burnings and the neutrino 
cooling. Furthermore, the thermal conductivity 
A n is dependent on the magnetic field. The MHD 
induction equation is supplemented as 



0. 



(5) 



From the mass conservation law ([2]), We obtain 
the mass flux m n as an invariant. 



pv = const. 



(6) 



On the other hand, equations governing the 
flow along the magnetic field are listed as follows. 



dx 



[pv 



P) 



0. 



(7) 



d (v 2 , 

pv dx-[Y + h 



dT 



where A p is the thermal conductivity parallel to 
the magnetic field. The contributions of magnetic 
field in equations (|7J| and © vanish because of the 
uniform magnetic field throughout the flame. 

The mass flux becomes an invariant again, but 
takes a different value. 



in. 



pv = const. 



(9) 



Because A p is not dependent on the magnetic 
field (see Section 2.3), all the above equations are 
independent of the magnetic field. 

We ignore the diffusion of ions. Therefore the 
Lewis number for each nuclei is assumed to be zero 
in the flames concerned here. 

2.2. Equation of state 

The equation of state used for this problem con- 
sists of ideal gas of ions, partially degenerate elec- 
trons together with electron positron pairs, and 
photons, which are in thermal equilibrium. We 
make a numerical table of all thermal quantities 
as functions of chemical abundance of an isotope 
i (or a nuclide i) Yi, the degree of electron degen- 
eracy ry, and the temperature T. 



2.3. Thermal conduction 



We follow the procedure of lPotekhin (jl999h for 



evaluations of thermal conductivities A p and A n . 
Suppose the magnetic field lies along the z-axis, 
the electron thermal conductivities introduced in 
the previous sections can be expressed as compo- 
nents of a 2nd-rank tensor such as A p = (X e ) zz 
and A„ = (Ae)^, because the collision time of 
electrons, which is dependent on the direction, be- 
comes a tensor Ty given by 

(10) 
(11) 

(12) 



Tzz = T(), 
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where u> is the cyclotron frequency of an electron, 
To is the non-magnetic electron collision time. The 
other components are zero. The tensorial conduc- 
tivity A in the non-quantizing magnetic field is 
written in terms of the following 2nd-rank tensors. 

A o = fc 2 T(/ 2 -/ 1 / - 1 / 1 ), (13) 

where k is the Boltzmann constant, and the ten- 
sors /„ are expressed as 



(In)i 
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(14) 

where p is the chemical potential of an electron, 
e the kinetic energy of an electron, H the Planck 
constant, m G the rest mass of an electron, and p is 
the momentum. The relation between e and p is 
(e + m c c 2 ) 2 = (m c c 2 ) 2 + (cp) 2 , where c is the speed 
of light. The effective mass m* of an electron has 
been introduced as m* = m c + e/c 2 . /(e) is the 
Fermi-Dirac distribution function, 



/(e) = 



1 



l + cxp(e-^)/fcT' 



(15) 



In the limit of extrem ely degenerate electrons 
( Urpin fc Yakovlevlll980l) . A e is approximated by 



n x ir 2 n e k 2 T 



(16) 



where n e is the electron number density and ep is 
the electron Fermi energy. The ratio of the ther- 
mal conductivities is expressed as 

^ ~ - . (17) 

A P 1 + ^ro) 2 1 ' 
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Therefore the existence of a magnetic field reduces 
A n and A n < A p . As a consequence, Equation 
(HJ) suggests that a magnetic field suppresses the 
transverse laminar flame velocity. 

2.4. Nuclear reaction 

We calculate a nuclear reaction network com- 
posed of the triple-a reaction, (0,7) reactions, 
(a,p)(p, 7 ) reactions, 12 C + 12 C, 12 C + 16 0, and 
16q + 16q ; with re i evant 25 isotopes from 1 H 
to 56 Ni (see Table [T]). This reaction network is 
smaller than the one used by TW92. The reac- 
tion rates and the energy g eneration rates S ar e 
referred to REACLIB (e.g. ICvburt et al1l201tt . 
Rate equations except for 4 He and 12 C are rep- 
resented in the following form, 



DYj 
Dt 



v— 

dx 

= ^i?4(p,T)^ - Rii(p,T)Yi 
3 

1 



1 + s jk 
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(18) 



where Yi is the abundance of a nuclide i with the 
mass number A; related to the mass fraction Xi 
hyYi = Xi/Ai, , R\l the decay rate of a nuclide i 
to a nuclide j, and R^ is the rate of the reaction 
between nuclide i and j to produce a nuclide k. 
The rate equations for 4 He and 12 C include con- 
tributions from the triple-a reaction in addition 
to the right hand side of equation (fT8|) . Addition- 
all y, we considered the electron screening effects 



all y, we considered 
bv iltoh et all (|l990t ) 



3. Boundary conditions 

We solve the MHD equations in each direction 
(eq. [2]-[3 for the transverse direction and eq. [7]- 
[5] for the longitudinal direction), with the rate 
equations (|T8|) satisfying the following boundary 
conditions. 

As boundary conditions in the upstream, we as- 
sume the state of an unburned WD material. The 
initial composition of the WD is assumed to be 
A( 12 C) = 0.5, A( 16 0) = 0.5. For a given den- 
sity, we set the initial temperature to the value at 
which the heating rate due to the nucle ar burning 
is equal to the neutrino cooling rate (jltoh et al 



19961 ). The other boundary conditions d/dx = 
is required at x — > 00. This is an eigenvalue prob- 
lem for the mass flux m that satisfies the boundary 
conditions at the two points. We seek the value of 
m to satisfy all the boundary conditions by inte- 
grating the MHD equations with various values of 
Vq, the velocity at x = in the wave rest frame. 
When a trial value of m is too large, the tem- 
perature increases faster with increasing x than 
the true solution, eventually diverges, and never 
reaches x —¥ 00. On the other hand, when m is 
too small, the temperature eventually decreases to 
zero before reaching x — > 00. Thus we can find the 
true value of m between values of these two trends, 
and the corresponding vo becomes the flame veloc- 
ity. 

We obtain solutions with the densities in 
the unburned media po ranging from 10 7 to 
10 9 g cm -3 . 

4. Results 

4.1. Flame propagation along the mag- 
netic field 

First, we present results for flames propagating 
along the magnetic field, which is equivalent to 
flames in non-magnetic media. We make a com- 
parison of the resultant flame velocities with those 
derived by TW92. 

When the initial density is po = 10 9 gcm -3 , 
the eigen value is found to be m p ~ 3.1 x 
10 15 g cm -2 s _1 , and the flame velocity is vq ~ 
31kms _1 . Figure Q] shows the resultant profiles 
of density and temperature. The result of nuclear 
abundances is presented in Figure The carbon 
burning leads to the production of elements up 
to sulfur, which is consistent with the result of 
TW92. 

In Figure [3l the flame velocities are compared 
with the result of TW92. The flame velocities 
of our results arc smaller than those of TW92. 
In TW92, a smaller reaction network (including 
fewer number of nuclei and reactions) resulted in 
smaller flame velocities. Since TW92 used a reac- 
tion network larger than our network, the smaller 
flame velocities of our results may be partly due 
to this difference. It is also claimed in TW92 that 
the velocity derived from the eigenvalue method is 
smaller than from the dynamical equations. The 
difference is prominent at low densities. Different 
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reaction rates and thermal conductivities in our 
calculations might be another factor to deviate the 
flame speeds from those of TW92. 

4.2. Flame propagation across the mag- 
netic field 

Figures @] and [5] show the solution for a 
flame propagating in the direction perpendicu- 
lar to the magnetic field with the initial density 
po = 10 9 gcm~ 3 and the initial magnetic field 
Bo = 1 12 G. The corresponding eigen value be- 
comes m n ~ 1.8 x 10 15 g cm -2 s" 1 , and the flame 
velocity vo ~ 18kms _1 . The scale of the burning 
region becomes shorter than that of the flame pre- 
sented in Figures Q] and [2j The velocity becomes 
slower due to the reduction of the thermal con- 
ductivity by the magnetic field. Both of the scale 
and the velocity decrease with increasing B . 

To see where the magnetic field affects the prop- 
agation of a flame normal to the magnetic field, 
Figure [6] shows the distribution of the ratio of con- 
ductivities A n /Ap. Here A n is the conductivity in 
the direction normal to the magnetic field and A p 
is evaluated by using the thermal state at each 
position and neglecting the magnetic effects. The 
velocity of the flame is suppressed by a factor of 
~ 2 due to the magnetic field (see Figs. (JlJ and 
(j4|) ) . This means the suppression of the conductiv- 
ity by a factor of ~ 4 from Equation (JXJ) . Therefore 
Figure [6] indicates that the magnetic field in the re- 
gion where carbon burning takes place determines 
the flame velocity. 

4.3. Critical magnetic fields 

We calculate the flame velocities v n normal to 
the magnetic fields for various strengths Bq of 
fields and densities p . Figure [7] shows the depen- 
dence of the flame velocities v n on the strength of 
magnetic field for a few different densities. Equa- 
tions (JJ) and (fT7f imply v n /v p ~ Bq 1 in the 
limit of large Bq. Figure [7] shows this trend 
for each density. If we define the critical field 
strength B CI as the strength of the magnetic field 
that suppresses the flame velocity v n by a fac- 
tor of 2, the critical magnetic fields are esti- 
mated as B CI ~ 10 10 ,10 n , andl0 12 G for den- 
sities po = 10 7 , 10 8 , and 10 9 gem -3 , respectively. 
Therefore even when the magnetic pressure does 
not dominate entire pressure, it is possible that 



the magnetic field suppresses the flame propaga- 
tion through the thermal conduction. 

Figure |8] shows the dependence of the flame 
width S = (Tb — T u ) / \dT / dx| max on the strength 
of magnetic field. Here, Tb and T u are tempera- 
tures in burned and unburned media. \dT / dx\ max 
is the maximum value of temperature gradients in 
the wave front. The flame width decreases when 
the strength of magnetic fields increases above the 
critical value, which is the same trend as v n shows. 
This implies that the burning timescale (~ S/v n ) 
is not affected by the existence of a magnetic field. 

5. Conclusions and discussion 

We study the propagation of the laminar flames 
in MWDs. We calculate the laminar flame ve- 
locities parallel and perpendicular to the mag- 
netic field as the eigensolutions of steady MHD 
equations taking into account nuclear reactions 
and thermal conduction by electrons. From the 
results, we estimate the critical magnetic fields 
B cl for given densities in the upstream region 
relevant to the interior of a C+O WD. Conse- 
quently, we find that magnetic fields of the order of 
B cr ~ 10 12 G significantly suppress the flame ve- 
locity at the center of the WD (the density thereof 
is p ~ 10 9 gcm~ 3 ) with the mass close to the 
Chandrasckhar limit. 

Though only the electron thermal conduction is 
concerned in our calculations, photons also trans- 
port the heat generated by nuclear reactions. The 
photon thermal conductivity increases with de- 
creasing densities and increasing temperatures. 
On the other hand, the electron conductivity is not 
sensitive to these quantities, becau se of the high 
degeneracy. IChamulak et al.1 (|2007| ) calculated the 
propagation of nuclear flames in non-magnetic me- 
dia with the thermal conduction including photons 
as well as electrons. They found that photons are 
more efficient than electrons for the heat transport 
in lower density media (p < 7 x 10 8 gcm~ 3 ). 

Recent observations have discovered many 
MWDs with strong surface magnetic fields up 
to 1 x 10 9 G. The mass of the MWD with 
the strongest magnetic field discovered so far is 
~ 0.9 Mq. The corresponding central density is 
estimated to be p c ~ 10 7 gcm -3 . If this MWD 
were to increase its mass close to the Chan- 
drasckhar limit, carbon would be ignited at the 
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center with the density of ~ 10 9 gcm -3 . Mean- 
while the strength of the magnetic field would 
increase scaling with the density as B ~ p 2 / 3 un- 
less the electrical resistance dissipated the mag- 
netic energy. Thus B will be enhanced by a fac- 
tor of ~ 20. According to lYoshida &: Eriguchi 



(|2006l ). there exist equilibrium configurations of 
non-rotating magnetized polytropic gas with the 
index N = 3 that has the ratio of the field 
strengths at the center to the pole as high as 
~45. Therefore the magnetic field is expected to 
become 20 x 45 x 10 9 ~ 10 12 G at the center, com- 
parable to the critical strength. In addition, a prc- 
supcrnova WD has exper ienced convection heated 
by 12 C + 12 C reactions (jHillebrandt fc Niemever 
2000T ). It will amplify the local magnetic field, 



which makes it more likely that the magnetic field 
exceeds the critical strength. 

(l2004h 



Ghezzi et al 



calculated the turbulent 
flame velocities propagating in the MWD. They 
considered effects of magnetic fields only on the 
growth rate of the Rayleigh- Taylor instabilities. 
The turbulent flame velocity was estimated by 
fractal scaling. Here the laminar flame veloc- 
ity corresponds to the maximum wave length of 
the Rayleigh- Taylor instability and the turbulent 
flame velocity the minimum wave length. The 
minimum wave length depends on the magnetic 
field through the growth rate of the instability, 
while the laminar flame velocity was assumed to 
be independent of the magnetic field. We find that 
the laminar flame velocity is also affected by mag- 
netic fields even weaker than those that can change 
the growth rate of the Rayleigh- Taylor instability. 
If this effect found in t his paper is taken i nto ac- 
count, the procedure of lGhezzi et al. I (l2004h would 
suggest more asymmetric explosions. In this ap- 
proximation, magnetic fields slow down the nu- 
clear burning. Therefore the MWD can expand to 
a larger extent before the flame reaches the sur- 
face. As a result, the nuclear burning occurs in 
a lower density environment, which leads to the 
production of a less amount of 56 Ni. However, 
there are some non-linear effects that may lead to 
significantly different results. For example, there 
must be a feed-back of a suppressed flame velocity 
perpendicular to the magnetic field to the thermal 
conduction along the field. The residual heat from 
the nuclear burning that cannot be conducted in 
the direction perpendicular to the field may be 



conducted along the field and change the flame 
velocity in this direction. Change of nuclear burn- 
ing may alter the expansion of the WD compared 
to a non-magnetic explosion. To investigate such 
effects, we need to perform multi-dimensional cal- 
culations, which is beyond the scope of the present 
paper. 
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Fig. 1. — Temperature and density profile 
of a flame parallel to the magnetic field, with 
po = 10 9 gcm~ 3 . The eigenvalue m p ~ 3.1 x 
10 15 g cm -2 s _1 , and the flame velocity vo ~ 
Slkms" 1 . 



This 2-column preprint was prepared with the AAS IATfrjX 
macros v5.2. 



Table 1: Nuclear isotopes contained in our reac- 
tion network. The top panncl shows rcactants and 
products of (a, 7) reactions. The bottom pannel 
shows those of (a, p) reactions. 



(a, 7) 


a, 12 C, iB 0, 2U Ne, 24 Mg, 2 *Si, 32 S, 
36 Ar, 40 Ca, 44 Ti, 48 Cr, 52 Fe, 56 Ni 


(a,p) 


p, lb N, i9 F, 23 Na, 27 Al, 31 P, 3b Cl, 
39 K, 43 Sc, 47 V, 51 Mn, 55 Co 
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Fig. 2. — Distributions of chemical abundance of 
a flame parallel to the magnetic field, with po = 
10 9 g cm~ 3 . 
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Fig. 3. — Laminar flame velocities for different 
densities in the unburned region. 
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Fig. 5. — Distributions of chemical abundance of 
a flame perpendicular to the magnetic field, with 
po = f0 9 gcm~ 3 , B = 10 12 G. 
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Fig. 4. — Temperature and density profile of a 
flame perpendicular to the magnetic field, with 
9™™-3 r„ — 10 12 G. The eigenvalue 



Po 



= 10 M gcm- d , B = 



m n ~ f .8 x fO 15 gem 
ity — 18 kms" 1 . 



and the flame veloc- 
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Fig. 6. — Ratio of conductivities A n /A p as a func- 
tion of position x for a flame propagating normal 
to the magnetic field. Here A n is the conductiv- 
ity in the direction normal to the magnetic field 
and A p is evaluated by using the thermal state at 
each position and neglecting the magnetic effects. 
pa = 10 9 gcm- 3 , B = 10 12 G. 
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Fig. 7. — Velocities of laminar flames prop- 
agating across the magnetic fields as func- 
tions of field strengths for densities of po = 
10 7 , 10 s , 10 9 gcm" 3 . 
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Fig. 8. — Widths of laminar flames prop- 
agating across the magnetic fields as func- 
tions of field strengths for densities of po = 
10 7 , 10 8 , 10 9 gcm- 3 . 
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